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Abstract. We consider a jet-disk symbiosis model to ex- 
plain Gamma Ray Bursts and their afterglows. It is pro- 
posed that GRBs are created inside a pre-existing jet from 
a neutron star which collapses to a black hole due to mas- 
sive accretion. In our model we assume that the initial 
energy due to this transition is all deposited in the jet by 
magnetic fields, using fully well explored concepts from 
jets and disks in active galactic nuclei and compact active 
stars in binary systems. The observed emission is then 
due to an ultrarelativistic shock wave propagating along 
the jet. We show that a good agreement between model 
predictions and observational data can be obtained for 
systems with accretion rates as high as in the Galactic jet 
source SS433. Specifically, we are able to reproduce the 
typical observed afterglow emission flux and its spectrum 
as a function of time. 
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1. Introduction 

Gamma Ray Bursts (GRBs) have been a mystery for al- 
most 30 years. Recently, thanks to the Italian-Dutch satel- 
lite BeppoSax, it has been possible to detect for the first 
time their X-ray, optical, and radio afterglows. Many arti- 
cles have been published to report the main character- 
istics of GRBs, (e.g., Guarnieri & al. 1997; Piro & al. 
1998; Frail & al. 1997; Metzger & al. 1997; Gorosabel & al. 
1998; Kulkarni & al. 1998). By now we know that GRBs 
are isotropically distributed over the sky, at least three of 
them are cosmological, they show a time variability in the 
7-emission of the order of milliseconds, and long complex 
bursts. Many, but not all of them have an X-ray afterglow, 
and it seems that for most of them a host galaxy can be 
found. 

Modern theoretical attempts to interpret the data are 
based on ideas by Meszaros and Rees (1993), Meszaros 



(1994), Panaitescu and Meszaros (1998), as well as 
Paczynski (1986), and Paczyhski and Rhoads (1993). In 
these models a relativistic shock is caused by a relativis- 
tic fireball in a pre-existing gas, such as the interstel- 
lar medium or a stellar wind, producing and accelerating 
electrons/positrons to very high energies, which produce 
the gamma-emission and the various afterglows observed. 
The low level of associated radiation at other wavelengths 
limits the baryonic load of the emitting regions to very low 
amounts, and constrains the scale of the emitting region 
to lengths much larger than a neutron star. 

One serious question is whether the overall energetics 
of the fireball — assumed to be isotropic — are reasonable 
(Sari and Piran 1997; Dar 1997) or actually exceed the 
level given by any conceivable model of neutron star merg- 
ers or other stellar collapses. Another question is whether 
a fireball model that uses external or internal shock waves 
can solve the baryonic mass load problem to explain the 
7-emission, starting with an initial energy of 10 51 erg in 
the spherical shell rest frame. 

Many authors suppose the validity of the fireball model 
and provide evidence for its agreement with the observa- 
tions (Waxman 1997a, Waxman 1997b, Vietri 1997), but 
the observational results for GRB971214, requiring an ini- 
tial energy E ~ 10 53 erg (if the emission is isotropic), pose 
a serious challenge to existing models. 

Here we propose an anisotropic model where the 7- 
ray emission and the afterglow is produced inside a pre- 
existing jet and calculate the temporal evolution of the 
corresponding flux. 

2. GRB jet model 
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In our model GRBs develop in a pre-existing jet. We con- 
sider a binary system formed by a neutron star and an 
O/B/WR companion in which the energy of the GRB is 
due to the accretion-induced collapse of the neutron star 
to a black hole. 

The high speed in the energy flow in AGN-jets is gener- 
ally believed to be initiated by strong magnetic field cou- 
puglieseQmpifr- phng to the accretion disk of the black hole (e.g. Falcke & 
Biermann 1995, Romanova & Lovelace 1997, Romanova 
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et al. 1997). We use the same mechanism here. We as- 
sume that in this transition a large amount of energy is 
anisotropically released as Poynting flux along the polar 
axis. One may think of this process as a violent and rapid 
twist, such as occasionally debated as a possible cause for 
supernova explosions (Kardashev 1970, Bisnovatyi-Kogan 
1970, Lc Blanc & Wilson 1970, and more recently Bier- 
mann 1993). This energy release naturally initiates an ul- 
trarelativistic shock wave in the pre-existing jet. The emis- 
sion microphysics are as in existing fireball models. 

To fix the jet parameters, we use the basic ideas of 
the jet-disk symbiosis model by Falcke & Biermann (1995, 
1999). In this model the direction of the magnetic field is 
mostly perpendicular to the axis of the jet and the val- 
ues of the total particle number n 3 (relativistic electrons 
+ thermal protons) and magnetic field B 3 are calculated 
using the equipartition between magnetic field energy den- 
sity and kinetic plasma energy density in the umperturbed 
jet. Both, B 3 and n 3 , are a function of the jet ejection rate 
Mj e t ~ 0.05M d isk ~ (10~ 5 M Q /yr)M_ 5j , the equipartition 
parameter e and the bulk velocity of the jet (3 3 -f 3 — 0.3«o.3, 
where j 3 is the bulk Lorentz factor of the jet prior to the 
burst. The values of the particle number density and the 
magnetic field in the unperturbed jet are then given as 

n 3 {z 3 )~ 80.34 (M- 6j %3)iry*r7j cm" 3 (1) 



B 3 (z 3 ) = 1.05 (M^v^ie-^OZI^) x 

z^ 1 - Gauss (2) 

where 7 m .2 is the minimum electron Lorentz factor 7 em = 
1007 m . 2 °f the- relativistic electrons (assumed to be in a 
power law distribution), 6-i 3 = #j/(0.1rad) is the opening 
angle of the jet, and z 3 = 10 17 zi7j cm is the distance along 
the jet. z is the redshift. 

The advantages of using the propagation of an ultra- 
relativistic shock wave inside a jet can be summarized in 
the following two points 

• The initiation of the shock by magnetic fields does not 

necessarily involve considerable matter. It implies a 
low amount of baryonic matter, of order 5 x 10~ 8 M Q , 
in the jet. 

• The initial amount of energy deposited in the jet is 

of the order E§\ = £?/(10 51 erg), and we can ob- 
tain an apparent isotropic energy as high as E app — 
10 53 3 E 5 i9zl } erg. This fits well the requirements of 
GRB971214. 

The evolution of the shock Lorentz factor 7 s h with dis- 
tance Zj along the jet axis can be obtained using the con- 
servation of the unperturbed jet gas energy in the shock 
rest frame, E/2 = it9 2 zfj 2 h m p c 2 n 3 , that is 

lsh (z 3 ) « H.48 (El[ 2 MZll 2 vl /2 )z- l7 f. (3) 



The characteristic time scale to see the emission across 
the entire region when the shock has reached a distance 
z } - along the axis of the jet is given by t^ ^ = z 3 /(2^ 2 h c). 
Substituting this into the formula for 7 s h, we get the time 
evolution of the ultrarelativistic shock front 

z i {t)^2Mxl^\El[ 2 MZli 2 vl / l)t\ /2 cm (4) 
%h (t) » 6.84 (E^MZ^)^, (5) 

where t 5 = i/(10 5 s). 

A change in the emission properties occurs when the 
opening angle of our jet passes the relation l/j s h = 0- 3 . 
Using the formula of the bulk Lorentz factor as a function 
of time, we can see at what time t* this relation holds: t* ~ 
2.20 x 10 4 (i?5iM_5jWo.3)#iij s. Hence, after about 6 hours 
this limit is reached. Prior to this, the observed emission 
is limited by the Lorentz boost to a conical section of the 
shock front of angle l/7 s h- For the following calculation 
of the afterglow emission several hours after the burst we 
can therefore consider 1 /j s h > 9 3 . 

In the flow behind a steady shock front, relativis- 
tic particles are usually accelerated and magnetic fields 
can be amplified. After the shock, Bj ~ B\\ and B' ± ~ 
47 s h-B±, relative to the jet axis. The resulting magnetic 
field in the shocked plasma will have a strength of roughly 

1 /2 

B w (167g h + 1) B± and the shock wave compresses the 
magnetic field component perpendicular to the jet axis. 
The jump conditions for the density particle number give 
n (sf) ~ 47 sh nj, (De Hoffmann and Teller 1950, Marscher 
and Gear 1985). We use here the approximation for an 
ultrarelativistic shock so that 712/ 'n\ = 4712 ~ ^Zsh, where 
1 and 2 are related to the zone before and after the shock. 
This allows a straight and simple limit to nonrelativistic 
shocks. In our model, this corresponds to 

n^(t) = 2.78 x 10W /4 Mlg%3 /4 ) x 

9zl 3 tZ 5/i cm- 3 (6) 

B,^\t) = 10.24 {E^M^T) x 

(e- 1/2 0Zl 3 Z^)h 3/4 Gauss (7) 

where (sf) shows the quantities in the shock frame and 
(ob) the corresponding values in the observer frame. 

In our jet we consider an electron power law distribu- 
tion N(j c )d-f c = C c -fZ p dj c where 7 om < 7 C < 7 max , with 
a cutoff at a constant minimum electron Lorentz factor 
7 cm ~ 100. In fact, if p-p collisions inject a population of 
electrons in the unperturbed jet, then one would expect 
that in the unperturbed jet the electron Lorentz factor j c 
goes from 7 cm ~ 100 ~ m^jm c up to some large value 
(Falcke & Biermann 1995). 
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The equation for C has been obtained considering 
that in the shock frame the powerlaw electron energy den- 
sity is taken as a fraction 8 < 1 of the pre-existing rela- 
tivistic electron energy density and using the value p=2 

Cf\t) = 1.91 x 10 4 £(£ 5 Y /2 M!g %3 3 / 2 ) x 

{ezl 3 i m .2)h V2 cm ~ 3 (8) 

If we were to relate to the proton energy density in the 
shock, 5 < 20 ~ m p /(m c 7 cm ). 

To calculate the boosting factor of the transition from 
the shock frame to observer frame, we assume the angle 
between the jet-axis and the line of sight to the observer 
is # bs < l/7sh some time during the 7-ray burst. The 
critical synchrotron frequency after the shock front in the 
observer frame is given by ^ m « TTzTsh7em^%> where we 
consider the minimum electron Lorentz factor evolves in 
phase space with the bulk Lorentz factor 7 om x 7 s h : 

^W-^^^ 2 ^- 1 ^^:^^)^ 372 Hz (9) 

We assume isotropic emission in the shock rest frame. 
Only when we transform the radiation emitted into the 
observer frame, the photons are concentrated in the for- 
ward direction, lying within a cone of half- angle 1 /7 s h- For 
the afterglow, however, the solid angle is determined by 
the actual opening angle of the jet. 

At the critical frequency, ^ m , in our model the syn- 
chrotron cooling time is less than the dynamical time, 
z j/(7shc). In fact, in the shock frame 

tW/tW « 7.88 x 10- 3 (Ml 5 >o. 3 )(^ 2 - lj 7 m 2 2 )t5 (10) 

so the synchrotron time scale is shorter than the dynami- 
cal time scale for our standard parameters. 

Considering the cooling evolution of our electron power 
law spectrum, to calculate the flux we use the formula 

rt° h \t) = mtl^ 6 ^ where (^)p°=2 - ^e to- 
tal emission power per unit volume per unit frequency for 
an electron power law distribution with an index p = 2, 
(e.g. Sect. 6.4 of Rybicki & Lightman, 1979). In our model 
this corresponds to a flux in the observer frame given by 

F(° h \t) = 4.00 x 10- ls cW(t)BV 2 W(t)v-V 2 z 2 (t) x 

x " e h*h x 7^2- ( n ) 

Here, the quantity x v represents the thickness of the 
radiating shock front, extending behind the shock front to 
a point where the local v m drops below v 

x v = 4.70 x ^{E^MZl^vyi) x 

( e 3 M 2 7 m 3 2 /4 )^V 4 1/2 cm (12) 



where v u = ^/(10 14 Hz). In both equations (11) and (12) 
v is in the shock frame. 

Substituting the values of the equations (4), (5), (7), 
(8) and (12) in (11), and with v in the observer frame 

Fl ob \t) = 7.45 x lO-^^MI^Vs ) x 

7m,2£ ) 2~8 2 5^ 5/4 ^r4 1 er S cm^s^Hz -1 (13) 

where D is the luminosity distance, D 2 s.5 = -D/(10 28 5 cm) 
corresponds to a redshift of about z ~ 1.5, using q Q = 1/2 
and H = 65 km s _1 Mpc _1 . 

The flux at frequencies v below the critical frequency, 
is m , can be obtained using the flux F Vtn (t) at frequencies 
v > I'm, calculated at the frequency v m , that is F u<Urri = 
F Uin (t) x ■ The corresponding flux is given by 

F^l{t) = (2.84 x 10- 28 )(2.61)-^(£; 5 3 1 /4 - /3/2 Mr 5 1 / 4 x 
^ / 3 4 )(e 1 / 2+ ^ 2 ^f7- 3 2 /2 - 5/3/2 )^- 8 2 5 x 

f l/4 + 3/3/2^ (1 + z)1+P ^ cm -2 s -l Hz -l (14) 

For optically thin emission, (3 is equal 1/3 and the time 
dependence is i 3 / 4 . This suggests a gentle optical rise. 

Another important check is to calculate the syn- 
chrotron self-absorption coefficient a v in our model (Sect. 
6.8 of Rybicki & Lightman, 1979). Integrating it over an 
isotropic distribution of particles and using the value ob- 
tained for the coefficient C c , we find a v {t). The optical 
depth along the path of a travelling ray is equal to the 
product of the absorption coefficient a v (t) times the width 
of the shell (equation 12), and it will be equal to unity at 
the time 

t T=l = QA2x S^HE-^M^v-^e-^OZf x 

7*>r 4 28/15 - (is) 

3. Discussion and conclusion 

We now check our model against the observed data. Here 
we assume a fixed opening jet angle and a minimum elec- 
tron Lorentz factor, and the equipartition parameters e 
and S are referenced to unity. 

For GRB970208, using the equation (9), at a time of 
topt = t/(7.6 x 10 4 s), we obtain a frequency f m (t op t) ~ 
(3.95 x 10 14 )(1 + zT'El[ 2 {e-V 2 ez\^l)t-V 2 Hz. 

Substituting the same numbers into the equations 
for the flux (13), we have F„(t opt ) « 3.09 x 10~ 28 
SiE^MZ^ ^ 3 4 ) 7m , 2 D 2 - 8 2 5 K P \ Cp/ 4 org ernes' 1 
Hz" 1 , where v opt = v/(3A x 10 14 Hz). 

Again using GRB970208 as reference, we obtain for 
the flux in the X-band at a frequency of 10 18 Hz 
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and a time 3.6 x 10 4 s: F v (t x ) « 2.67 x lCT 31 ^^ 

For the optical depth in the optical band we obtain 
r v {t opt ) « 2.00 x 10- 13 5(£- 5/8 Ml 5 5 f % 15/8 )(e-V^ : ^ 

7m,2) v opt 2 topt^ ■> hence, within our model the optical 
emission region is optically thin, except very early. 

In our model, the condition l/7 s h < Oj corre- 
sponds to the early phase, when v m is actually in the 
gamma-ray regime. In this phase, the corresponding max- 
imum flux in the observer frame is -F m ax ~ 8.66 x 
lO^SE^D^Ozlj-ymM" 1 er S cm~ 2 s _1 . This is derived 
from using all available energy in energetic electrons, and 
redistributing it into emission; photons are generated by 
various emission processes, including pion decay, upscat- 
tered by inverse Compton emission, and redistributed in 
photon energy by pair opacity (Rachen & Meszaros 1998). 
The integration of this flux in time from 10~ 3 s to 10 s 
gives a value of about W~ 5 5E 5 iD2$ 5 0Zij"fm,2 erg cm~ 2 . 
This we identify with the initial gamma-ray emission. 

Concerning the variability at times when l/j s h > 6 } , 
we note that in observed jets, one often finds inhomo- 
geneities on the scale of a few jet-diameters parallel to 
the axis, and down to some fraction of the diameter per- 
pendicular to the axis. However, any variability derived 
from such inhomogeneities is smeared out by arrival time 
differences for the observer. Therefore, using the jet diam- 
eter as a reference scale, this smearing limits any temporal 
variability 

^>(7^) 2 - 0A7(El?M:l!\tt)t^ 2 (16) 

Inhomogeneities on transverse scales smaller than the 
jet diameter can shorten this. This variability may explain 
the complex features of the optical rise in GRB970508 
(Galama et al. 1998). 

The way in which our model can reproduce the ob- 
served properties of GRBs depends strictly on two sets of 
parameters, one set which is verifiable, because it derives 
from known active binary stars, and a second set which 
characterizes the explosion. The explosion depends on the 
initial shock energy £51, and the equipartition parame- 
ter S. On the other hand, the binary set is composed of 
the mass loss rate M_ 5 j of the jet, the jet opening angle 
9j, the bulk velocity in the jet «o.3 an d the equipartition 
parameter e. To produce shock waves and the observed 
emission in the jet, it is necessary that the initial amount 
of matter in the shock is comparable with the mass in the 
jet. This implies that the binary system necessary in our 
model has to have a super Eddington accretion rate to 
get the required mass loss rate. In our Galaxy, one such 
binary system is known, SS433 (Murata and Shibazaki, 
1996), for which M disk ~ (2 x 10" 4 M Q /yr). In our model, 
we have used the parameters characteristic of this system, 
implying that there are system like SS433 where the cen- 
tral object indeed is a neutron star. 



In the context of the Falcke & Biermann jet-disk 
model, SS433 is a radio-weak jet-disk system, while here 
we have used their radio-loud model; but it has been 
noted that systems such as GRS1915+105 and perhaps 
also SS433 can switch from radio-weak to radio-loud in 
the terminology of Falcke & Biermann (1995, 1999). Con- 
sidering that the jet ejection rate is bound to the disk ac- 
cretion rate by the relation Mj et ~ 0.05Mdi s k, we find that 
for an initial bulk Lorentz factor r o = 10 , the mass in 
the shock M s h « 5 x 10 _8 Mq is equal to the mass present 
in the jet. Therefore, SS433 might be a good candidate to 
explode as the next violent GRB in our Galaxy. 

To summarize, our model can explain the initial 
gamma ray burst, the spectrum and temporal behaviour 
of the afterglows, the low baryon load, an optical rise, and 
do all this with a modest energy budget. Moreover, this 
GRB model is developed within an existing framework for 
galactic jet sources, using a set of well determined param- 
eters. 

Obviously, we have simplified in many places, using 
a naive version of shock acceleration, only strong shocks, 
a conical jet geometry, etc., but the good agreement with 
the data we obtain within our framework shows that more 
detailed calculations may be worthwhile. 
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